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It has become common practice to model large spin ensembles as an effective pseudospin with total 
angular momentum J — N x j , where j is the spin per particle. Such approaches (at least implicitly) 
restrict the quantum state of the ensemble to the so-called symmetric Hilbert space. Here, we argue 
that symmetric states are not generally well-preserved under the type of decoherence typical of 
experiments involving large clouds of atoms or ions. In particular, symmetric states are rapidly 
degraded under models of decoherence that act identically but locally on the different members 
of the ensemble. Using an approach [Phys. Rev. A 78, 052101 (2008)] that is not limited to the 
symmetric Hilbert space, we explore potential pitfalls in the design and interpretation of experiments 
on spin-squeezing and collective atomic phenomena when the properties of the symmetric states are 
extended to systems where they do not apply. 
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I. INTRODUCTION 

For a variety of fundamental and technological reasons, 
there is considerable interest in studying quantum fluc- 
tuations in the angular momentum of large atomic/ionic 
spin ensembles [H[2l[3l|4l|5j[6llTJ|Sl[9l [10l [11]. From a 
theoretical perspective, modeling such systems is compli- 
cated by the fundamental property of quantum mechan- 
ics that the Hilbert space ,34?n describing TV spin-j par- 
ticles grows exponentially with the number of particles, 
dimJ^v = (2j+l) Jv . As aresult of exponential scaling, it 
has become common practice to look for dynamical sym- 
metries that reduce the effective dimension of the spin 
ensemble by restricting its state to a manageable sub- 
Hilbert space [TSJ [T3] . One then makes inferences about 
the properties of the large ensemble based on those of the 
sub- Hilbert space. But, of course, the validity of such in- 
ferences depends critically on how well the actual spin 
system respects the symmetries used to formulate the 
reduced-dimensional description of its quantum state. 

Although limited exceptions exist [T2J [TJ] , most work 
to date on reducing the effective dimension of large spin 
systems has focussed on the symmetric group [T51[TS1HB] : 
the sub-Hilbert space ^ C -^n spanned by iV-body 
states that are invariant under the permutation of par- 
ticles — IV'), \ip) <= ^S- In theory, the sym- 
metric group provides a model of experiments that 
cannot distinguish between particles during any por- 
tion of state preparation, manipulation or measurement. 
For spin-1/2 ensembles, the dimension of the symmet- 
ric group grows only linearly in the number of spin- 
1/2 particles, dim^| = N + 1 <C 2 , making it ex- 
tremely amenable to simulation and analysis. Yet, the 
symmetric states still exhibit interesting multi-particle 
phenomena, such as entanglement [13J, spin-squeezing 



|17j and zero-temperature quantum phase transitions 
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This favorable trade-off between manageable size and 
non-classical behavior has made the symmetric group 
the sub-Hilbert space of choice for analyzing large spin 
ensembles — indeed, any approach that models a large 
spin ensemble as a collective pseudospin of size J = N x j 
[3 El ESI [161 [22J [23J EU E3 [261 E7] is grounded at least 
implicitly in the theoretical underpinnings of particle ex- 
change symmetry [IS]. To justify using the symmetric 
Hilbert space as a realistic model, two key assumptions 
are generally made: 

• Assumption 1: The degree of spin polarization 
that is achieved in practice (such as by optical 
pumping and possibly additional purification) is 
sufficient to prepare the ensemble into a state that 
is well-described by a nearly-pure symmetric state, 
and ideally by a spin coherent state. 

• Assumption 2: Symmetric states are nearly pre- 
served under low to moderate levels of decoherence, 
at least of the variety typically encountered in prac- 
tice, such as that due to spontaneous emission of a 
far-detuned probe laser. 

Furthermore, it is generally taken to be true that rea- 
sonable laboratory efforts to achieve homogeneous cou- 
pling to the electromagnetic fields used to manipulate 
and measure the ensemble correspond to conditions well- 
approximated by permutation invariance. Under these 
assumptions, several key properties of symmetric collec- 
tive states, reviewed in Section (II A 2 1, have played a 
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central role in the design and interpretation of experi- 
ments involving large spin ensembles: 

• Interpretation 1: Spin-polarized states, such as 
those obtained by optical pumping, exhibit min- 
imum uncertainty in angular momentum observ- 
ables transv erse to the direction of polarization 
(AJj_ 4 ) = \/N/2 with respect to the Heisenberg- 
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Robertson inequality (in units where h = 1) 

(AJ ±1 )(AJ ±2 )>\([J ±1 ,J ±2 })\/2. (1) 

• Interpretation 2: Classical noise, or the uncer- 
tainty that results from a classical mixture of spin 
eigenstates, grows faster than y/~N and linearly in 
TV for the worst case. Projection noise scaling that 
grows faster than y/N can be used to diagnose the 
presence of classical uncertainty in the ensemble. 

These properties are such fundamental characteristics of 
the symmetric states that (at least some) research groups 
have been known to train their members to view a linear 
increase in spin polarization coinciding with a square- 
root increase of spin-projection noise [2] with atom num- 
ber as a laboratory signature of a spin coherent state. 
These misconceptions have very likely led to the mis- 
characterization of spin-squeezing in all but perhaps the 
most recent experiments on spin-noise reduction in large 
atomic ensembles [251 |2"9"1 IBUl EH] . 

A. Symmetric versus Collective Decoherence 



operators. It has thus become common practice [22J [25] 
to study decoherence in spin ensembles by approximating 
Eq. Q by its associated collective process 



£ c {S]p = 



SpS^ - i (&Sp + p&S 
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Eq. ([5]) is more amenable to analysis and simulation be- 
cause it preserves the (N + l)-dimensional symmetric 
states. But, it is not always a good physical model. In 
atomic systems, for example, a typical source of deco- 
herence comes from spontaneous emission, yet collective 
radiative processes only occur under stringent conditions 
such as superradiance from highly confined atoms [32] 
and some cavity-QED or spin-grating settings [I]. Even 
in these cases, the extent to which N atoms behave as a 
single point-particle dipolc moment is imperfect at best. 
Under typical experimental conditions, where an atomic 
or ionic ensemble is coupled to a free-space laser probe 
and the average interatomic spacing is not small com- 
pared to the laser wavelength, Eq. (|5| is just simply the 
inappropriate model of decoherence [331 151] . regardless 
of how much easier it is to analyze. 



For the symmetric Hilbert space Jtfg to remain an ac- 
curate description of a spin ensemble's state (provided 
that the initial state is an clement of J#s), the system's 
dynamics must be generated by completely symmetric 
collective processes: processes that are themselves per- 
mutation invariant and thus expressible in terms of col- 
lective operators 



(2) 



Such operators apply the same single-particle operator 
s G 0ti(2) to each atom in the ensemble, where sW = 
li ® • • • <8> l„_i ® s n ® l„+i ■ • • ® • • • In acts non-trivially 
only on the n th particle. As such, Eq. Q is explicitly 
permutation invariant by construction. 

Unfortunately, many of the decoherence models most 
appropriate for large spin ensembles cannot be described 
as collective symmetric processes even when the decoher- 
ence acts identically on each particle. Consider, the open 
system dynamics governed by the master equation 



dp(t) 
dt 



-jC[s]p{t), 



(3) 



where decoherence acts with the same rate 7 but locally 
on every member of the ensemble via the Lindbladian 
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The (TV + l)-dimensional symmetric-group Hilbert space 
J^s is not preserved under such dynamics, as the Linblad 
superoperator cannot be expressed in terms of collective 



B. Main Results 

In this paper, we argue that the statistics of collec- 
tive angular momentum operators in large spin ensembles 
are not well-predicted by the behavior of the symmetric 
group. Many of our results fall in stark contrast with the 
conventional wisdom surrounding large spin ensembles: 

• The uncertainty in collective spin obervables for the 
completely depolarized state of N spin- 1/2 parti- 
cles scales as \/N with the number of particles, a 
scaling that is analytically equivalent to that of a 
pure spin coherent state. 

• For large ensembles, optical pumping does not pro- 
duce an approximately pure symmetric state even 
at high levels of spin polarization. For example, 
even with an optical- pumping efficiency of 99.9%, 
the purity of an ensemble with N ~ 10 6 spin- 
1/2 particles is vanishingly small, tr[p 2 ] ~ 10 -409 , 
while its overlap with the symmetric group is about 
10 -205 . Both decrease exponentially with TV. 

• For partially polarized ensembles (e.g., incomplete 
optical pumping) , the uncertainty in transverse col- 
lective spin observables (AJj_i) scales as \fN with 
the number of particles, while the polarization ( jy) 
scales linearly in N. Thus, essentially every state 
of the ensemble corresponding to incomplete opti- 
cal pumping exhibits the same scaling behavior as 
an actual spin coherent state. 

• Even if one could prepare a pure initial coher- 
ent state, even small levels of decoherence rapidly 
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transform the ensemble state into one that is ex- 
tremely mixed and very poorly described by a sym- 
metric state. In fact, we predict that an ensemble 
with N ~ 10 5 particles which has decohered by 
20% (its polarization has dropped to 80% that of 
the initial coherent state) has a rather small purity, 
approximately 1Q -92 ' 630 . 



A. Generalized Collective States 

In the "irrep basis," arbitrary pure states of the spin 
ensemble are expressed as 



^ cj t M,i\J,M,i), 



(14) 



J.M.i 



II. COLLECTIVE STATES OF THE ENSEMBLE 

Consider an ensemble of N identical spin- 1/2 particles 
described by the single-particle Pauli operators cn n ) = 

(a^\ ay n \ ai^ ) and corresponding angular momentum 

:■(») _ ~(») 



operators ja"' = fra /2 [38]. The joint Hilbert space 
for the entire spin ensemble JT = Jf W <g> • • • <g> W 
has dimension dim(.Jf) = 2^, and arbitrary pure states 
of the ensemble can be expressed in the tensor product 
basis 



^ ] c m 1 , 



\mi,m 2 ,...,m N ) 



(6) 



which still requires 2 N coefficients [refer to Eq. ( 13 )] . Of 
course, simply transforming to the irrep basis does not 
change the effective dimension of the Hilbert space, but 
it suggests the symmetry that was used to develop the 
concept of generalized collective states in Ref. |12j . Such 
states are described by the sub- Hilbert space Jtf?c C J#n 
spanned by TV-particle states that are indistinguishable 
across the d J N degenerate irreps for each total angular 
momentum J. This generalized permutation symmetry 
cj,M,i — cj,M,v j Vi, i' makes it unnecessary to distinguish 
basis kets \J,M,i) with respect to their irrep label. By 
defining effective basis kets \J,M) on each total- J irrep 
block, the generalized collective states are 



where the basis states \mi, . . . , tun) = \\-,rn\)\ ® ••• (E 
l!,mjv)jv are simultaneous eigenkets of [j^™* 1 ] 2 and ji™' 1 : 

m 2 \ 



mi, 
i (n) K, 



,m N ) = j n (j n + l)|mi, . . 
,m N ) = m n |mi,...,mjv) 



,m N ) 



(7) 
(8) 



Each particle in the ensemble transforms separately 
under rotation such that = (R)]® N , where 
S>i(R) is the spin-1/2 rotation operator parameterized 
by the Eulcr angles R — (a,/3,7). Expressed in the 
tensor-product basis, the [&^(R)]^ N provide a reducible 
representation for the rotation group but can be decom- 
posed into irreducible components (irreps) 



9{R)= 



J=J„ 



0i^(i?) 



via the total spin eigenstates 

J 2 \J,M,i) = J(J+l)\J,M,i) 
j z \J,M,i) = M\J,M,i) 



(9) 



(10) 
(11) 
and 



with the collective spin operators J q = \ S n =i ^ 
J = mod(iV/2, 2), . . . , N/2. For each total angular mo- 
mentum J, the quantum number i = 1 
guishes between the 



, d J N distin- 



iV!(2J+ 1) 



(N/2 - jy.(N/2 + j + iy. 



(12) 



degenerate irreps with total angular momentum J |35j . 
It is readily shown that the degeneracy function satisfies 



N/2 

£(2J+1)4 

,7=0 



(13) 



He) = J2 cj > m \ j ' m }> 

J,M 



(15) 



with the rcscaled coefficients 



/l 

CJM = A -7J- ^ Cj,M,i, 



(16) 



where the summation is over the d^j copies of the irrep 
with total angular momentum J . Under this symmetry, 
dim(Jff c ) = (N + 2) 2 /4 (for N even) scales only quadrati- 
cally with the number of particles in the ensemble. While 
not as convenient as the linear dimensional scaling of the 
symmetric group, the 0(N 2 ) scaling of the generalized 
collective states is still a vast improvement over expo- 
nential scaling and is sufficient to allow simulations with 
at least a hundred or so particles. 

When studying decoherence and other open-system 
dynamics of a spin ensemble, it is necessary to work with 
the density operator of the system, rather than a state 
vector. The collective state density operator is defined as 
the direct sum over the reduced density operators pj for 
each total- J irrep block [T^] 



PC 



= ©^ = EE Pj,M;J,m'\^M){J,M'\. (17) 



J M,M 



As defined, the collective density operator restricts 
against coherence between irrep blocks. It is shown in 
Ref. [T5] and in Section ( III I that symmetric maps exhibit 



a type of super-selection property, which prevents them 
from generating coherences between irrep blocks. Gen- 
eralized collective states are therefore sufficient to model 
any dynamics of the form in Eq. Q, provided that the 



initial state satisfies Eq. (17 1 
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1. Irrep Populations and Purity 

The structure of generalized collective states can be 
analyzed by considering the fraction of the population 
that resides within each total- J irrep block 



PJ = tr[pj] 



(18) 



where pj is the reduced density matrix defined in Eq. 
(17 1. The overlap of a generalized collective state 



with the symmetric group is therefore given by — 
tr[/5/v/2]- Another important distinction between sym- 
metric and generalized collective states is that the collec- 
tive states can be mixed over total- J irrep blocks even 
if all of the reduced density operators pj are internally 
pure. From Eq. (16 1 the purity of the full density opera- 
tor p is given by 

Purity = tr[p 2 ] = £ jytr^]. (19) 



2. The Symmetric States 

The symmetric states previously considered for large 
spin ensembles are a special case of the generalized col- 
lective states [12]: J%s C J#c is spanned by the maximal 
angular momentum manifold: cj m = for J ^ N/2 and 

d N/2 

N 



1. Since the symmetric states only have sup- 
port on the maximum- J irrep, their trace vanishes on all 
irrep blocks except for J = N/2, tr[pj] — Sj_ N / 2 , thus 
providing a simple test to determine whether a collective 
state is also symmetric. The symmetric states can be de- 
fined equivalently as the manifold of states that can be 
reached from the maximum j z eigenstate \N/2, N/2) «-> 
1 1 , . . . , 2) using only maps generated by collective oper- 
ators, Eq. ([2]) (the z-polarized state \N/2,N/2) is clearly 
permutation invariant). 

The spin coherent states are a special case of the 
symmetric collective states, defined by the manifold of 
states that are simply connected to the z-polarized state 
\N/2,N/2) by a rotation 



\0,cf>) = ®(0,<f>)\N/2 > N/2). 



(20) 



It is well-known, and readily shown, that the expecta- 
tion value along the direction of spin polarization Jq^ = 
9(0, (p)J z &(0, 0) is given by 



(0,^^0,0) = (N/2,N/2\J z \N/2,N/2) 



N 



(21) 



with (AJg^) — while the transverse expectation values 
vanish for the spin coherent state 

(J ±1 ) = l -(N/2,N/2\(J + + J_)\N/2,N/2)=0 (22) 
(J ±2 ) = l -(N/2,N/2\(J + -J_)\N/2,N/2)=0, (23) 



but the variances do not 
(A 2 JO - ^(N/2,N/2\J+J_\N/2,N/2)-0 = ~. (24) 



The scaling of the spin projection noise (AJj_.) = yN/2 
for a coherent state is the basis for Interpretation 1, 
described in Sec. HI 

The completely mixed state of a spin system with total 
angular momentum j is given by p = i2j+i/(2j + 1), and 
therefore the completely mixed symmetric state is 



remixed 



1 



iAT/2©6 A r /2 _ie---e6, (25) 



which has the property that it is completely depolarized 
with respect to all collective spin operators 

1 



-tr[J„] = 0. 



N + 1 

The variance in collective spin observables 



N + l 



12 



(26) 



(27) 



is a direct consequence of the permutation-invariancc 
constraint tr[p7\r/2] = 1- This linear scaling of the spin 
projection noise for the mixed state /5^ ixcd is the basis 
for Interpretation 2, described in Sec.|T] 



3. The Completely Mixed Collective State 

When permutation-invariance is lifted (retaining in- 
variance only over the degenerate copies of irreps), the 
completely mixed state of the N spins generalizes to 

1 d ' N h 1 

Pmixed = S ^TTT = ( 28 ) 

J a N i=1 LJ + 1 Z j 

That is, for each irrep contribution to the direct sum, the 
elements of the density operator are given by the ratio 
of the degeneracy of that irrep to the total dimension of 
the HUbert space, 



Pj,m-j,m 



i — 2 dff5M,M' 



(29) 



iV!(2J+l)<5 Mi M< 



2 N (N/2 - J)\(N/2 + J + 1)! ' 



precisely as would be expected. Normalization of the 
completely depolarized state is readily verified using Eq. 
(13 1. Once again, the expectation values of all collective 
angular momentum operators vanish 



(Ja) = 0, 

but their variance does not 



(30) 



<A 2 ia) = ^tr 



N/2 

E 

J=0 



3 • 2 



iY 



N 



(31) 
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In fact, the uncertainty of all collective spin observables 
with respect to the completely mixed state is quantita- 
tively identical to that of the spin coherent state 



(AJ a ) = 



III. SYMMETRIC DYNAMICS AND 
COLLECTIVE-STATE PRESERVING 
PROCESSES 

The super-operator 

£f>C = E E PJ-M;J,M'fiiM' 



(32) 



./ M,W 

will preserve collective states L : Mq — > if its action 

f J MM ,=C\J,M){J,M>\ (34) 

can be expressed in the irrep basis in such a way that it 
does not distinguish between degenerate irreps f^M 1 = 
/mm" Operators that transform simply with respect to 
the rotation group 

§ = 0& = E E sJ m,m>\J,M)(J,M'\, (35) 

J J MM 1 



including the collective angular momentum operators J a 
and all collective operators 5 = X)n=i f° rme d from 
s( n ' S su(2), satisfy the requirement of invariance over 
degenerate irreps by construction. 



But as discussed in Section (IA), processes that are 



only symmetric over local single-particle super-operators, 



£ s [«]p=^£(«)[«(")]p = ^§W^("))t 



(36) 



do not transform simply under rotations. Our present 
work is made possible by results from our previous 
demonstration that any symmetric local map of the form 



in Eq. ( 36 1 can be brought into the form of Eq. ( 33 ) and 



therefore preserves collective states [H] . For the su(2) 
operator 



s ■ a = sq! + s+(T-| 



(37) 



expressed in the basis {<r_ , er+, <t z , 1}, the action of Eq. 
( 36 1 can be constructed as 



f( IM ,=s-g(J,M,M')-^ 
from the tensor operator 



9qr 



J2^ n) \J:M)(J,M'\ai n 



)t 



(38) 



(39) 



The elements of g can be derived recursively [T^] to give 



aj,m aJ.M' 

q A r 

g qr = 2J x 



(40) 



a J N +1 (2J+l) 
d J N (J + l) 



\J,M + q)(J,M' + r\ 



b jm b jm> J 

2Jd? - \ J -^ M + l){J-l,M' + r\ 

a J N +1 D J > M Di> M ' 

q - fj \J +l,M + q) (J + 1 , M + r\ 



2(J+ l)d 



N 



where the reduced degeneracies are given by a ,] N 



(33) ^2j'=j d J N and the coefficients are defined as 



A J ± M = +V(JT M)(J±M+ 1) 



J,M 



= M 



J,M 



= ±y/(JTM)(JTM-l) 



B J Z ' M = \/(J + M)(J — M) 



D 



J,M 



tV(J ± m + ± m + 2 ) 

Df M = y/(J + M+l)(J-M+l). 



N — 

(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



The three terms in Eq. ( 40 1 arise from two types of pro- 



cesses: (Term 1) transitions that occur between M levels 
within a single J irrep; and (Terms 2-3) transitions that 
couple neighboring irreps with A J = ±1. It is this cou- 
pling between irreps that prevents maps of the form in 
Eq. Q from preserving symmetric states and that makes 
collective models of decoherence inadequate for modeling 
spin ensembles under most laboratory conditions. 



IV. EXAMPLES 

We have found simulations of large spin systems to be 
an invaluable tool for studying the properties of symmet- 
ric decoherence. Even though the effective dimension of 
the generalized collective states grows faster that that 
of the symmetric group, 0(N 2 ) rather than O(N), it is 
still possible to run simulations over a sufficient range to 
make both qualitative and quantitative predictions. As 
such, we have performed simulations aimed at addressing 
the following specific questions: 

1 . "Is it possible to prepare a large spin system into a 
state that is well-approximated by a spin coher- 
ent state, and thus a symmetric state?" Time- 
evolving the density operator for a spin system un- 
der a model of optical pumping enables us to ana- 
lyze the purity and irrep structure of the system as 
it is spin-polarized from an initial mixed state, in- 
cluding what happens for incomplete polarization. 

2. "Do symmetric states remain a good model of spin 
ensembles subject to limited decoherence?" Time- 
evolving the iV-particle density operator under a 
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symmetric model of spin depolarization enables us 
to study the relationship between the expectation 
value and uncertainty of collective angular momen- 
tum operators as well as the irrep structure of the 
state as it decoheres from an initial coherent state. 

3. "Is the practice of approximating symmetric de- 
coherence models with their associated collective 
processes justified if only expectation values and 
uncertainties of collective operators are of inter- 
est?" Time-evolving the system under an entan- 
gling Hamiltonian and contrasting the effect of the 
different decoherence models allows us to compare 
their collective statistics. 



A. Partial Polarization of the Spin Ensemble 

To determine whether symmetric states, and in par- 
ticular spin coherent states, provide a good description 
of a large spin ensemble subject to optical pumping, we 
considered the symmetric polarizing channel 



dp(t) 
dt 



lC s [ 3+ ]p{t), 



(47) 



which describes the effective spin-1/2 dynamics that arise 
when radiative excited states are adiabatically eliminated 
from atoms with two ground states under conditions 
where the atoms are coupled to a circularly-polarized 
laser field [36 . It is readily shown that the steady state 
corresponding to the symmetric polarizing channel is the 
spin-coherent state \0 — 0, <fr — 0), i.e., the state that is 
polarized along the positive z-axis, with (J z ) = N/2. 

Under typical laboratory conditions, however, optical 
pumping does not achieve complete polarization; pump- 
ing falls short of reaching the steady state of Eq. (47) 
[3TJ|37]. Figures [l|al) and[TJbl) plot the time-evolution 
of the collective expectation value ( J z ) and uncertainties, 
(AJ X ), (AJ y ) and (AJ Z ), as the spin ensemble evolves 
from a completely mixed initial state, Eq. ( 28 I for N — 50 
and N = 100 particles. As expected, the spin polariza- 
tion (J z ) increases monotonically from its initial value 
of zero, coinciding with a decrease in (J z ) as the sys- 
tem progresses toward the maximum- J z eigenstate. The 
transverse uncertainties, (J x ) and (J y ), are constants of 
the motion: beginning at \J N /2 for the completely mixed 
state and remaining at y/ N/2 at all times as the system 
progresses toward the |0,0) coherent state. 

Figures ^&2) and[l}b2) show the irrep decomposition 
and purity of the spin ensemble as it is gradually po- 
larized under the dynamics of Eq. ( |47| ). The reduced 
traces pj — tr[/5j] is shown for each of the total- J ir- 
rep blocks, for J = 0,1,..., N/2 (for clarity, only the 
irrep blocks with J close to J ma x — N/2 are labeled on 
the plot). The initial completely mixed state has an ex- 
tremely small overlap with the symmetric group and a 



purity that is exponentially small in N. As the dynam- 
ics proceed, population is gradually transferred to irreps 
with increasing angular momentum. Furthermore, the 
progression of the state to higher total- J irrep blocks is 
apparently slower for N — 100 particles than for N = 50 
particles. For N — 50, the maximum- J irrep (symmetric 
group) begins to show a non-negligible population when 
the spin polarization is approximately 80% of N/2. For 
N = 100 particles, however, the symmetric group does 
not begin to be populated until nearly 90% spin polar- 
ization. The behavior of the purity [dashed lines in Figs. 
\^a2) and[l|b2)] is more dramatic. Even at 98% spin 
polarization, the state of the ensemble is far from pure: 
tr[/5 2 ] < 0.4 for N = 50 and tr^ 2 ] < 0.2 for N = 100. 

The apparent decrease in purity and overlap with the 
symmetric group for a given level of spin polarization as 
the number of particles is increased is explored further 
in Fig. [3] For each value of A, the state of the system is 
evolved under Eq. ( 47 I from a mixed state at t = until 



the time when the fractional spin polarization 



/ 



N/2 



(48) 



achieves a target value. The corresponding state is then 
analyzed to determine its overlap with the symmetric 
group pn/2 = tr[/3jv/2] and its purity tr[p 2 ]. Fig. p[a) 
plots the results for fractional polarizations / = 92%, 
95%, and 98% over the range 4 < N < 120. As can be 
seen from Fig. [3ja) , both the overlaps with the symmetric 
group (solid lines) and the purities (dotted lines) decrease 
exponentially with the number of spins ./V over the range 
of N that could be analyzed. Given the consistency of the 
simulation data as a function of N, it seems reasonable 
to extrapolate the results to higher values of N by fitting 
the data to an exponential form: 



tr[/5 2 



10 



and tr[/3jv/ 2 ] « 10 



° P N 



(49) 



Values of the exponents r)° p and rf/v for various frac- 
tional polaizations / are listed in Table [j] The results 
are quite dramatic, suggesting that even at very high 
levels of spin polarization, such as / = 99.9%, the pu- 
rity and symmetric overlap achieved by optical pumping 
in typical experiments are both vanishingly small, eg., 
tr[p N / 2 ) ~3x 10~ 21 for N = 10 5 . Thus, it seems reason- 
able to conclude that symmetric states are vastly inade- 
quate for describing such ensembles. 

Figure |2jb) shows the scaling of the mean polarization 
(J z ) and the transverse uncertainties, (AJ X ) and (AJ y ), 
as a function of A for different levels of optical pump- 
ing efficiency. As expected, the mean polarization scales 
linearly with N. For incomplete optical pumping, its 
value is degraded with respect to the maximum value by 
the factor /. That is, (J z ) = /A/2. More surprisingly, 
perhaps, is that the transverse uncertainties, (A j x ) and 
(AJ y ), are always equal to \/N/2 regardless of the degree 
of spin polarization. This result illustrates that there is a 
fundamental flaw in the laboratory practice of identifying 
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FIG. 1: (color online) Simulation of a model of spin-polarization dynamics, corresponding for example to optical pumping 
computed for N — 50 (plots al and a2) and N = 100 (plots bl-b2) spin- 1/2 particles. Beginning from the completely-mixed 
state Eq. (281, the spin ensemble evolves under the symmetric z-axis polarizing channel given by Eq. (47 1. As the dynamics 



proceed (al and bl), the mean polarization (J z ) increases while its uncertainty (AJ Z ) decays. The transverse uncertainties, 
(AJ X ) = (AJ y ) = N/2 are a constant of the motion. The evolution of the J-irrep block traces (plots a2 and b2) clearly 
shows that despite ( J z ) quickly approaching maximum polarization, the J < J max irrep blocks are still highly populated and 
the state is quite mixed. 



a spin coherent state simply from scaling behavior: spin 
polarization that scales as N coinciding with transverse 
uncertainty that scales as \/~N . Rather, such an iden- 
tification is only possible provided with a high-quality, 
independent measurement of N . 



B. Decoherence from a Spin Coherent State 

To assess whether symmetric states provide a good 
model of large spin systems subject to decoherence, we 
considered the dynamics 



TABLE I: Fitted values of the scaling exponents in Eq. (491 
for the purity tr[p 2 ] and overlap with the symmetric group 
tr[/3jv/2] that is partially polarized to a fixed level of spin- 
polarization / beginning from an initial completely-mixed 
state. 



V° P P 



80.0% 
90.0% 
95.0% 
98.0% 
99.0% 
99.9% 



0.0856200 
0.0424700 
0.0209700 
0.0084510 
0.0041670 
0.0004082 



0.0454300 
0.0218100 
0.0106200 
0.0042460 
0.0020890 
0.0002042 



dm 

dt 



(50) 



of an initial spin coherent state subject to the symmetric 
depolarizing channel 



C s - 

'-DP — 



(C s [j x } + C s [j y }+C s [j z ] 



(51) 



As discussed in the introduction, the symmetric depolar- 
izing channel acts identically but locally on each spin in 
the ensemble, which is in contrast to the collective analog 
of Eq. (51 1, and given by and collective 



c s DP = cF [X] + cF [j y ] + cF [j,] ) ■ (52) 
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FIG. 2: (color online) Demonstration that the symmetric group provides an extremely poor model of a large spin ensemble 
that is incompletely polarized under the symmetric polarization channel Eq. (47 1 beginning from an initial completely mixed 
state. Plot (a) shows the scaling of the purity tr[p 2 ] and overlap with the symmetric group tr[/5jv/2] for systems that have been 
polarized to 92%, 95% and 98% of (J™ ax ) = N/2. Both the purity and overlap with the symmetric group are seen to decrese 
exponentially in N. Plot (b) shows that the spin-polarization ( J z ) grows linearly in N while the transverse uncertainties grow 
as %/N even for systems that are incompletely polarized. 



Figure [3] plots the time evolution of the expectation 
value (,/ z ) and uncertainties, (AJ X ), (AJ y ) and (AJ X ), of 
the collective spin operators for ensembles consisting of 
N = 50 and N = 100 particles beginning from the ini- 
tial z-polarized spin coherent state \6 = 0, </> = 0). As 
expected, the expectation value (J z ) decreases in time 
while the uncertainty (AJ Z ) increases; however, the un- 
certainties (AJ y ) and (Aj x ) are constants of the motion, 
in contrast to the behavior that would be observed under 



Eq. (52) 



The irrep block traces are plotted as a function of time 
in the bottom panels of Fig. [3] for N = 50 and N = 100 
particles. It is evident from the plots that the symmetric 
group quickly becomes a poor description of the state 
of the spin ensemble: the trace tr[p^/ 2 ] °f the maxi- 
mum J-irrep block quickly decays. For N = 50 particles 
at 95% polarization tr[p2s] < 0.6, indicating that much 
of the population has been transferred to lower J-irrep 
blocks. For N — 100 at 95% polarization, the departure 
from the symmetric group is even more dramatic, with 
tr[/?5o] < 0.35. As N becomes larger, this behavior be- 
comes more pronounced and even minimal decoherence 
produces significant deviation from the symmetric states. 



Analogous to the fitting procedure described in Sec. IV A 



the purity and symmetric group overlap can be extrap- 
olated to higher numbers of particles according to the 
exponential fits 



tr[p 2 



10" 



and tx[p N / 2 ] w 10" 



(53) 



Values of the exponents rf^ and r]^ p for various deco- 
herence levels (measured by the remaining fractional po- 
larization /) are listed in Table [II] A gain, the results 



are dramatic, suggesting that even at low levels of de- 
coherence, the remaining purity and symmetric overlap 
become exponentially small. Thus, it also seems reason- 
able to conclude that symmetric states are vastly inade- 
quate for describing a spin ensemble subject to even small 
amounts of decoherence even if it were possible to prepare 
an initial coherent state. The scalings of the symmetric 
overlap and purity are plotted as a function of N in Fig. 
[5] for various levels of depolarization. 

Figure [4] depicts the structure of the density operator, 
expressed in irrep-block basis, for various levels of deco- 
herence. For the sake of clarity, the figure was generated 
for a rather small number of atoms N — 16, however, 
we have verified that the qualitative results generalize 
to higher N. As the dynamics proceed from the initial 
spin coherent state toward the depolarized state, irrep 
blocks with lower total angular momentum, J < J max 
become populated. The steadily decreasing value of the 
spin polarization ( J z ) is therefore a result of two mecha- 
nisms: deocherence within each irrep block, and mixing 
between the blocks. Throughout this process, the trans- 
verse uncertainties of the collective spin operators, (AJ X ) 
and (Aj y ) do not increase, even though the uncertainties 
in individual irrep blocks do. 



1. Projection Noise of the Depolarized State 



We compared the scaling of the uncertainty in the col- 
lective spin observables J x , J y and J z , for the steady 



state solutions to the symmetric Eq. (51 1 versus collec 



tive Eq. ( 52 1 mode models of decoherence as a function 
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FIG . 3: (color online) As the state of a spin ensemble decoheres from a coherent state under the symmetric depolarizing channel 
Eq. (51 1, the mean spin-polarization ( j z ) decreases from its maximum value of N/2; however, the transverse uncertainties, (A j x ) 
and (AJ y ) remain constant with the value y/ N/2. As the ensemble decoheres, population is transferred out of the symmetric 
group and into total-J irreps with J < J max , as indicated by the evolving irrep traces tr[pj(i)], as the purity of the state tr[p 2 ] 
(dotted line) decreases accordingly. Comparison of the decoherence dynamics for N = 50 (a) versus N = 100 (b) particles 
suggests that the state of the system leaves the symmetric group more rapidly as N increases. 



TABLE II: Fitted values of the scaling exponents in Eq. (531 
for the purity tr[p 2 ] and overlap with the symmetric group 
tr[/3jv/a] that has decohered to a fixed level of spin-polarization 



/ beginning 


from an initial spin-coherent state. 




/ 


Vp p 


nf 


80.0% 


0.09263 


0.04944 


85.0% 


0.07182 


0.03768 


90.0% 


0.04946 


0.02552 


92.0% 


0.04165 


0.02138 


95.0% 


0.02552 


0.01296 


98.0% 


0.01719 


0.008685 


99.0% 


0.008685 


0.004365 



of the number of particles N. To do so, we determined 
the steady state density operator p ss corresponding to 
the decoherece dynamics by solving 



0. 



(54) 



In practice, this is accomplished by expressing both the 
quantum state p and the superoperator L in their associ- 



ated Liouville representations, where p is an 0(N 2 x 1)- 
dimcnsional column vector and L is an 0(N 2 x iV 2 )- 
dimensional sparse matrix. The steady state density op- 
erator is then given by the eigenvector associated with 
the A = eigenvalue of C. 

Our results are illustrated by Fig. [6^ a) . The uncertain- 
ties of the collective spin observables (A J a ) for the steady 
state of the symmetric depolarizing decoherence superop- 
erator Eq. (51) as a function of the number of particles 
N are shown by the circles. As expected, for a com- 
pletely depolarized steady state, all of the collective spin 
uncertainties are identical (A j x ) = (AJy) = (AJ Z ). Fur- 
thermore, the uncertainties scale as yN with the number 
of particles, verified by the agreement of the data points 
with the predicted uncertainty scaling of the completely 
depolarized state derived in Eq. (32). For comparison, 



the uncertainties of the collective spin observables (AJ a ) 
were also computed for the steady state solutions to the 
collective depolarizing superoperator Eq. (52 1 as a func- 
tion of the number of particles [squares in Fig. [6][a)]. 
Again, all uncertainties are equal for the completely de- 
polarized steady state; however, the scaling with N is 
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FIG. 4: (color online) A fully polarized state undergoing no decoherence is confined to the symmetric space corresponding to 
the highest J-irrep block. As the state becomes decohered under the depolarization channel Eq. < |51[ ), population from the 
highest J-irrep is transferred to lower J-irrep blocks. Simulations show that even for N = 16 particles, as seen above, mild 
amounts of decoherence ~ 10% significantly deplete the J = 8 irrep block and the state is driven far from the manifold of 
symmetric states. At 30% decoherence, only a vestigial trace population remains in the J = 8 irrep block. 



linear, in accordance with Eq. (27 1 
Figure Ifjj 



shows the computational resources re- 
quired to simulate generalized collective states. The prac- 
tical limitation to the maximum value of N that could 
be analyzed was determined by the required to compute 
the A = eigenvector of Cfj P . Two hours were required 
to do so for N — 120 even though the memory required 



to store the Liouville representation of Eq. (51 1, which 
was only on the order of 30 MB for N = 120 (although 
the swap-space consumed by the eigensolver was at least 
5 GB for the N = 120 calculation). 



C. Decoherence and Dynamical Spin Squeezing 

As a final example, we compare and contrast the dy- 
namics of a spin system that is subject to Hamiltonian 
evolution generated by a collective operators as it un- 
dergoes symmetric versus collective models of decoher- 
ence. Specifically, we compare evolution under the mas- 
ter equation 



dp(t) 
dt 



-i[H,p(t)] 



l£ S D pp{t) 



(55) 



versus the master equation 
dp(t) 



dt 



= -i[H,p{t))+iC c DP p{t) 



for the "counter-twisting" Hamiltonian [17] 
H = -iX ( J? - J 2 



(56) 



(57) 



that has previously been used to study spin squeezing 
within the symmetric group. 

Dynamics were simulated for a variety of values of N, 
A and 7 beginning from a z-polarized spin coherent state 
\6 = 0,(f) = 0). Figure [7j a) plots the time-evolution 
of the collective expectation value (J z ) and transverse 
uncertainties, (AJ X ) and (AJ y ), for the specific case of 
N — 50, A = 1/50 and 7 = 4/50. As expected, the mean 
polarization ( J z ) decreases over time in all three cases. In 
the absence of decoherence this apparent depolarization 
is a byproduct of the increased uncertainty in J x , i.e., 
the anti-squeezing, even though the state remains pure. 
When decoherence is added, increased depolarization is 
observed, as expected [comparison of the dotted, dashed, 
and solid lines for (J z ) in Fig. [7J a)]. As can be seen, this 
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depolarization is noticeably more pronounced under the 
model of collective decoherence than for the symmetric 
model. The time evolution of the transverse uncertain- 
ties, (A J x ) and (A J y ), also shows that the symmetric 
model of decoherence degrades the collective statistics 
less so than does collective decoherence. 

Each taken alone, neither the uncertainty reduction 
(A j y ) nor the depolarization ( J z ) assesses the utility of 
the spin system for precision measurement [3 [TO] . Figure 
|7jh) plots the time-evolution of the squeezing parameter 



(Jx 



(58) 



which provides a metric for characterizing the sensitiv- 
ity of spin-resonance measurements relative to that of a 
coherent state (for which £ 2 = 1). It is evident from the 
plot that, even in the presence of decoherence, the squeez- 
ing parameter can drop below £ 2 = 1. Even though it 
would be technically incorrect to refer to the state of the 
ensemble as a "spin-squeezed state" (as such a state as 
typically defined constitutes a relatively pure symmetric 
state) [17] , it would appear that a metrological improve- 
ment over an actual coherent state is possible even in the 



presence of symmetric decoherence under Eq. (51 1 



V. CONCLUSION 

We have identified a number of flaws inherent in using 
the qualitative properties of symmetric states for mod- 
eling the behavior of large spin ensembles. Even in the 
most state-of-the-art laboratory settings, experiments in- 
volve atom numbers ranging from 5 x 10 3 — 10 7 . Under 
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these conditions, it is uncommon to find optical pumping 
efficiencies better than 95%-98%, and in many cases the 
degree of spin polarization may be much worse. How- 
ever, even in the best examples of spin polarization, we 
have found that the resulting state is not well described 
by a pure spin coherent state even if the expectation val- 



120 



FIG. 5: (color online) The purity tr[p 2 ] (dotted lines) and 
symmetric group overlap trfpjv/a] (solid lines) are plotted for 
a system that has decohered by 2%, 5% and 15% of (J™ ax ) = 
N/2 from an initial spin coherent state under Eq. (|51[). 
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FIG. 6: (color online) Plot (a) compares the uncertainty in 
the collective spin observables (AJ q ) for the steady states 
of the symmetric depolarizing channel Eq. (51 1 versus the 



collective depolarizing channel Eq. (52 1 as a function of the 
number of spin-1/2 particles, TV. In the case of collective 
decoherence, the calculated steady state of £§p agrees with 
Eq. ( |25[ ) and the calculated uncertainty (squares) matches 
that of the completely-mixed symmetric state, with (A J a ) = 
\J N(N + 2)/12 (dotted line). For the symmetric decoher ence 
model C% P , the calculated steady state agrees with Eq. (281 
and the calculated uncertainty (circles) matches that of the 
completely-mixed generalized collective state, with (AJ a ) = 
\fN /2 (solid line). In the latter case, the uncertainty scales 
identically to that of a spin coherent state, despite that it 
is completely mixed. Plot (b) indicates the computational 
resources required to find the steady state as a function of N: 
the memory required to store the Liouville superoperator for 
Cup (squares, right-side axis) and the time required to find 
the A = eigenstate p ss (circles, left-side axis). 
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FIG. 7: (color online) Simulations of N = 50 particles evolving under the counter-twisting Hamiltonian H = — i\( J+ — j£ ) with 
either collective, symmetric, or no decoherence. Plot (a) shows that although the expectation ( J z ) and uncertainty (A j x ) evolve 
in a qualitatively similar manner for the three models, the J y -uncertainty does not: the symmetric and no decoherence models 
exhibit a decrease in (AJ y ), but the collective model increases. As a result, the symmetric and no decoherence models exhibit 
spin-squeezing while the collective model does not, as evidenced in the plot of the squeezing parameter £ 2 = N{A 2 J y }/( J z ) 2 in 
plot (b). 



ues of collective spin operators achieve values that are 
approximately those of a coherent state. Furthermore, 
even when sufficient care is taken in the laboratory to re- 
duce decoherence to minimal levels, highly mixed states 
with little overlap in the symmetric group are inevitably 
produced. As a result of these findings, we conclude that 
greater care must be exercised when interpreting exper- 
iments on large spin systems using scaling laws inferred 
from the properties of symmetric states. 
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